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Using numerical renormalization group techniques, we study static and dynamic properties of a family of 
single-channel Kondo impurity models with axial magnetic anisotropy DS^ terms; such models are appropriate 
to describe magnetic impurity atoms adsorbed on non-magnetic surfaces, which may exhibit surface Kondo 
effect. We show that for positive anisotropy D and for any spin S, the systems behave at low temperatures as 
regular Fermi liquids with fiilly compensated impurity spin. The approach to the stable fixed point depends 
on the value of the spin S and on the ratio D/t'^\ where T^'' is the Kondo temperature in the absence of 
the anisotropy. For S = 1, the screening occurs in two stages if D < T^'; the second Kondo temperature is 
exponentially reduced in this case. More generally, there is an effective spin- 1/2 Kondo effect for any integer 
S* if -D < rj^?' and for any half-integer S if D > T^^ . For negative anisotropy D, the system is a non-Fermi 
liquid with residual anisotropic exchange interaction. However, the presence of transverse magnetic anisotropy 
E{S'^ — Sy) restores Fermi-liquid behavior in real systems. 

PACS numbers: 75.30.Gw, 72.10.Fk, 72.15.Qm, 68.37.Ef 



I. INTRODUCTION 

As the sizes of the active components in magnetic devices 
shrink, the ratio of the surface area over bulk volume steadily 
increases and the importance of surface and interface effects 
grows. In the past, the focus has mostly been on the mag- 
netic properties of clean surfaces, thin films, and interfaces 
in multilayers ''^••'•"', but now magnetism of nanostructures is 
emerging as a new subfield with the aim of exploring and 
modifying the exchange interaction at the nanometer level. 
Given that the symmetry in the surface area is always reduced 
as compared to the bulk, more attention needs to be given 
to possible additional surface-induced magnetic anisotropy 
effect a'i^'^'^-^-^'^''"''' . Recently, magnetic anisotropy on sur- 
faces studied at the atomic level has become a new focus of in- 
terest with potential applications in ultra-dense magnetic data 
storage. Using scanning tunneling spectroscopy, the magnetic 
anisotropy can be measured on the single-atom level by an- 
alyzing the splittings of magnetic excitation peaks as a func- 
tion of the applied magnetic field^iiS. Magnetic anisotropy 
constants were found to be very large, often in the mcV/atom 
range^'"'-'^. 

Single magnetic impurity atoms adsorbed on surfaces or 
buried in the near-surface regions can be described using 
quantum impurity models'-'-'"' with anisotropy terms'"-'^. The 
dominant anisotropy terms are of the form D{S-rif', where n 
is a vector along a privileged axis that is commonly chosen to 
be the z-axis'^ '^. The most important contributions seem to 
arise from the local spin-orbit coupling, as the reduced sym- 
metry leads to unequal hybridization of the different rf-orbitals 
of the impurity atom.'^''^. For an impurity buried in the near- 
surface region, D depends on the distance from the surface 
and oscillates in sign'**. 

Numerical calculations and experiments show that the mag- 
netic anisotropy strongly depends on the microscopic details. 
Depending on the impurity species and on the substrate, dif- 



ferent directions of the privileged axis may be found. Density 
functional theory calculation predict, for example, that per- 
pendicular anisotropy is common for single adatoms on metal- 
lic surfaces'"^. 

Spin-excitation spectroscopy using a scanning tunneling 
microscope (STM) has revealed a rich variety of different 
classes of behavior. For example, for Fe (spin-2) adatom on a 
thin CuN decoupling layer on a Cu surface, the anisotropy is 
negative with the z-axis in the surface plane, along the di- 
rection of nitrogen atoms, while for Mn (spin-5/2) on the 
same surface, the anisotropy is weaker but with an out-of- 
plane z-axis, however D is still negativsiS. No anomalies at 
the Fermi level (surface Kondo effect^SiSii^^S^"*) are observed 
for D < 0'°. For Co (spin-3/2) on CuN layer, the anisotropy 
D is positive and signatures of the Kondo effect are observed 
in tunneling spectra^^. These results strongly motivate a study 
of the entire family of the Kondo impurity models for various 
values of spin and for different signs and magnitudes of the 
magnetic anisotropy. Due to the low symmetry at the surface, 
it appears sufficient as a first approximation to consider cou- 
pling to a single conduction channel (the one with the largest 
exchange coupling constant); it is unlikely that coupling con- 
stants to two orthogonal channels would be equal under these 
conditions. 

Similar impurity models with large anisotropies also ap- 
pear in the context of single molecular magnets, e.g. Mni2 
or Feg. It was possible to attach metallic contacts to these 
molecules allowing for electron transport measurement a^^-^^ . 
The molecules exhibit very large spin values (S = 10 for 
Mni2) and a negative-Z? axial anisotropy term DS^, so that 
ground states with two polarized large-spin states are favored 
in the absence of quantum spin tunneling due to the transverse 
anisotropy E{Sl - S^) tsrm^^^2^. 

Since the ultimate miniaturization limit appears to be stor- 
ing information in the spin states of single magnetic atoms, 
it is important to study decoherence by possible magnetic 
screening effects through nearby conduction electrons (Kondo 
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effect), since the goal is to have long-lived large-spin magnetic 
impurity states. The appropriate model is thus the impurity 
Kondo model with additional anisotropy terms that we intro- 
duce in Secin] In Sec.|III]we relate these models by means of 
scaling equations to the Kondo models with XXZ anisotropic 
exchange coupling to the conduction band. Numerical renor- 
malization group results are presented in Sec. |IV| where we 
compare calculations for S — 1, S — 3/2 and S — 2. 
S = 3/2 is a typical representative for all half-integer-spin 
models, 5 = 2 a typical representative for all integer-spin 
models, while S = 1 is a special case with some additional 
features when compared to higher integer spins. We present 
results both for static properties (magnetic susceptibility, im- 
purity entropy) and for dynamic properties (spectral functions, 
dynamic magnetic susceptibility). This section also contains 
a discussion of the various fixed points, focusing on the non- 
Fermi-Uquid behavior in the D < models. 



II. MODEL 

A general approach to study the anisotropy effects of ad- 
sorbed impurity atoms would start with a multi-orbital An- 
derson model for the d-orbitals with a suitable hybridization 
function for the coupling with the substrate and with electron- 
electron repulsion and Hund coupling parameters extracted, 
for example, from density functional theory calculations. The 
next step would then consist of obtaining an effective Kondo- 
like model using a Schrieffer- Wolff transformation. We aim, 
however, to explore only the effects due to local magnetic 
anisotropy effects. We thus study a very simplified Kondo 
impurity model defined by the following Hamiltonian: 
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Here operators c^a describe the conduction band electrons 
with momentum k and spin a S {Tj i}- In numerical cal- 
culations, we will for simplicity assume a flat band efc = Wk, 
where the dimensionless momentum (energy) k ranges from 
— 1 to 1 and W is the half -bandwidth; the density of states is 
thus constant, p — 1/{2W). Furthermore, the effective ex- 
change constant J is defined as 
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where Jk.k' describe exchange scattering of the conduction 
band electrons on the impurity, and the spin-density of the 
conduction band electrons s is defined as 
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If Jk.k' = J is constant, then s can be interpreted as the spin- 
density of conduction band electrons at the position of the im- 
purity. For an impurity adsorbed on a surface, this clearly 
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Figure 1 : (Color online) Energy level diagrams for the anisotropic 
magnetic impurities with spin 1, 3/2 and 2 as a function of the abso- 
lute value of the anisotropy, |D|. Wide (blue) lines represent pairs of 
states, narrow (red) lines represent non-degenerate states. 



cannot be a good approximation as the exchange interaction 
Jk.k' is inevitably anisotropic in the reciprocal space; never- 
theless, since we focus mostly on the role of the local mag- 
netic anisotropy terms, we will not pursue this issue in this 
work. We furthermore assumed that the exchange interaction 
is isotropic in the spin space. 

In the anisotropy term of the Hamiltonian, 



DSt + EiSl 



(4) 



the first term is the axial magnetic anisotropy, while the sec- 
ond is named the transverse magnetic anisotropy. The axis 
z is by convention chosen such that is maximized, while 
axes xy are oriented so that £' > 0. If the axial anisotropy 
term is negative, < 0, we call the z axis an easy axis; if, 
however, _D > the anisotropy is said to be hard-axis (or pla- 
nar). Axial anisotropy term D leads to the splitting of mag- 
netic levels, Fig.[T] In the case of easy-axis anisotropy, Z) < 0, 
the moment tends to be maximal in size with two degener- 
ate states pointing in the opposite directions, i.e. Sz — ±5. 
Since the conduction band electrons can only change the im- 
purity spin by one unit during a spin-flip scattering event, the 
exchange scattering rate is expected to be strongly reduced 
for large enough \D\. In the case of hard-axis anisotropy, 
> 0, the single level = is favored for integer 
S and the doublet Sz — ±1/2 for half-integer S; in both 
cases, the impurity moment will be compensated at low tem- 
peratures, since the doublet can be screened via the conven- 
tional spin-1/2 Kondo screening. It should be noted that 
E{Sl - Si) = E/2 [{S+f + [S-f] term, i.e. the trans- 
verse anisotropy, mixes levels with different values of Sz- 
This term, even when E is small, plays an essential role in 
the easy-axis case^; in the hard-axis case, however, it only 
leads to a small correction. 

The isotropic Kondo model (i.e. the D = E — 
limit) has been intensively studied by a variety of tech- 
niques, such as numerical renormalization group22ii2i2ii^, 
Bethe Ansatg^^i3A38-39,40,4i^ ^jj^j ^jj^g^ methodsiSiMi. These 

studies have uncovered the physics of underscreening: a sin- 
gle conduction channel can screen at most 1/2 unit of spin, so 
that at low temperatures the impurity remains magnetic with 
a residual spin of S* — 1/2'*''. It was found that for S > I the 
approach to the stable Fermi-liquid fixed point is slow (loga- 
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rithmic); this behavior was characterized as that of a singular 
Fermi liquid^"*'^^ . 

The anisotropic Kondo model with D ^ (and with ad- 
ditional anisotropy in the exchange coupling) was studied by 
the Bethe-Ansatz technique'*^''*^ . It was shown that anisotropy 
can induce a quantum critical point (i.e. non-Fermi-liquid be- 
havior). The applicability of the Bethe-Ansatz approach is, 
however, limited to a set of models with restrained parameters 
thus a direct comparison with the model considered in this 
work is not possible. Nevertheless, there is qualitative agree- 
ment in that we also find a non-Fermi-liquid state for D < 0. 

The anisotropic Kondo model was also previously studied 
by the numerical renormalization group technique in Refs.[28i 
and I29I with the focus on high values of spin S and easy-axis 
anisotropy as appropriate for molecular magnets (we note that 
the definition of D in the cited works differs in sign from 
ours). Where comparisons can be made, our results agree with 
theirs. 

Finally, we mention that related models and physical ef- 
fects are also studied in the context of transport spectroscopy 
of quantum dots and impurity clusters. The high conduc- 
tance in the case of underscreening in spin-1 quantum dots 
and the two-stage Kondo screening were discussed in Refs.i48l 
and Furthermore, the two-stage Kondo screening is also 
found in the case of a single conduction channel in multiple- 
quantum-dot structures, in particular in the case of side- 
coupled quantum dots^"'^'-^^-^-'. High-spin states can be ob- 
tained in systems of multiple impurities if the conduction- 
band-mediated exchange interaction is ferromagnetio^i^, 
which may occur for small inter-impurity separations. Two- 
stage Kondo screening and the closely related "singlet-triplet" 
Kondo effec t^^'^^'^^'^^'^°'^^ have both been intensively studied 
experimentally^^iSi^i^^i^. 



III. SCALING ANALYSIS 

There is a close relationSSi^ between the Kondo model with 
the DS^ magnetic anisotropy term and the Kondo model with 
the XXZ anisotropic exchange constants and j ^69,70,7i.72,73 
defined by the Hamiltonian 



Hk = JzSzSz + J± [SxSx + SySy) 



(5) 



To be specific, we consider here the 5 = 1 anisotropic Kondo 
model with both types of the anisotropy (XXZ exchange and 
the DS^ term, but E — 0). Taking into account that the en- 
ergy of the 5^ = ±1 intermediate states is higher by D as 
compared with the energy of the Sz = state, the following 
scaling equations are obtained^liS^iZS 
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Here the scaling parameter is I — — In W; it runs from to 
positive infinity as the energy scale is reduced. Note that here 



we denote the half-bandwidth by W (this departs from the 
NRG convention of denoting it by D). We have introduced 
dimensionless running coupling constants by absorbing the 
density of states p: j± = pJ±, jz = pJz, while d is mea- 
sured in the units of W: d — D/W. These scaling equations 
differ from those in Ref. |6^ only by the d-dependent factors 
in Eqs. (|6|l and (|7]), which arise due to the energy shift by D 
in the 6*2 = ±1 states. 

From Eq. (O it follows that d will rapidly grow in absolute 
value*^, since the d term on the right hand side is first order, 
while the j'^ and terms are second order, hence smaller. 
The importance of the axial anisotropy is thus always rein- 
forced at lower energy scales. Furthermore, starting from ini- 
tially equal bare coupling constants, jz = j±,a non-zero d in- 
duces XXZ exchange anisotropy by unequal renormalization 
of J 2 and j±, as seen from Eqs. (|6]l and (|7]l. At low energy 
scales, the physical properties are qualitatively the same irre- 
spective of the origin of the anisotropy. From Eq. ([8]) we may 
also anticipate that the D > case corresponds to a > 
exchange anisotropic Kondo model, while the D < Q case to 
a J I < J2 model^. 



IV. NUMERICAL RESULTS 

Calculations were performed using the numerical renor- 
malization group techniqu e ' 3, 1432,33,74,7s y^jji^jj consists of a 
logarithmic discretization of the continuum of states of the 
conduction band electrons, followed by mapping to a one- 
dimensional chain Hamiltonian with exponentially decreasing 
hopping constants. This Hamiltonian is then diagonalized it- 
eratively by taking into account one new lattice site at each 
iteration. We have used discretization parameter A = 2.5 
(A = 2 for spectral function calculations) and truncation with 
the energy cutoff at i^cutofi = lOci^Ar, where uj^ is the char- 
acteristic energy scale at the A^-th iteration step. As a fur- 
ther precaution, truncation is always performed in a "gap" of 
width at least O.OIcjat, so as not to introduce systematic er- 
rors. To prevent spurious polarization of the residual impurity 
spin in the D < Q case due to floating-point round-off er- 
rors, it is helpful to symmetrize the energies of states which 
should be exactly degenerate (i.e. the ±5*2 pairs). In all calcu- 
lations, we have taken explicitly into account the conservation 
of charge Q; in calculations with only the DS^ anisotropy 
term (E = 0), a further conserved quantum number is Sz- 

In all calculations we set pJ ~ 0.1, therefore the Kondo 
temperature in the isotropic case, is the same for all val- 



ues of spin 
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Extracting T^^ directly from the NRG results we obtain a 

more accurate value of = 1.16 x lO^'^W. (We use Wil- 
son's definition of the Kondo temperature^i2i2Z.) 
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A. Static properties 

We first consider the static (thermodynamic) properties, in 
particular the impurity contribution to the magnetic suscep- 
tibihty, Ximp, and the impurity contribution to the entropy, 
5inip. The first quantity, when multiplied by the temperature, 
is a measure of the effective magnetic moment at a given tem- 
perature scale, while the second provides information about 
the effective impurity degrees of freedom. 

In Fig. ID we plot the thermodynamic properties of the 
anisotropic Kondo model with the axial DS^ term (E — 0) for 
three different values of the impurity spin S and for a range of 
the anisotropy strength. As long as the temperature is larger 
than the anisotropy D, the curves follow closely the results 
for the isotropic case (black curves). At lower temperatures, 
the behavior of the system strongly depends on the sign of 
D. For negative D (easy-axis case), the impurity spin remains 
partially unscreened at low temperatures; for small \ D\, the ef- 
fective moment is slightly above the isotropic strong-coupling 
fixed-point value with spin S -1/2 [i.e. fi ^ {S^ - l/4)/3] 
and it saturates with increasing ji^l at the value of /i = S*^. 
The low-temperature impurity entropy is In 2 for all negative 
D. This suggests that the impurity behaves as a residual two- 
level system with states |+) and |— ), which in the large- jZ?! 
limit become equal to the \Sz = +S) and \Sz = —S) states. 
The scattering of the low-energy electrons on this residual im- 
purity states is discussed in App. [A] For spin S* > 1, the 
effective moment always increases as the temperature T is 
decreased below \D\, which reflects the progressive freezing 
out of the magnetic levels other than the maximally aligned 
Sz = States. For spin 1, this is still the case for large 
enough \D\, but for lower \D\ the temperature dependence 
becomes monotonic. The boundary between the two differ- 
ent behaviors can be approximately located at jZ?! ^ T^k' ■ 
The reason appears to be that for \D\ < the 5* = 1 

Kondo effect had already significantly screened the impurity 
spin yielding a residual spin 1/2 by the time the anisotropy 
begins to be felt; the only effect of the anisotropy is then to 
prevent the screening process from completing, which gives a 
residual spin value slightly above 1/2. 

For positive D (hard-axis case), the impurity is non- 
magnetic at low temperatures for any value of the impurity 
spin S. Unlike in the easy-axis case, we see a larger variety of 
possible behaviors in the approach to the fixed point depend- 
ing on the spin and on the anisotropy strength. We discuss 
S ~ 3/2 and S = 2 first; these are characteristic representa- 
tives for all half-integer and integer spin cases. 

For half-integer spin and small D < the curves just 

follow the results for the isotropic case until T ^ D, when 
spin states other than the Sz = state of the residual (inte- 
ger!) — 1/2 impurity spin freeze out and the system ap- 
proaches the non-magnetic ground state exponentially fast. 
For sufficiently large D > the high \Sz\ states of the 

original impurity spin freeze out exponentially fast atT ^ D, 
this time yielding a Sz = ±1/2 doublet which then under- 
goes spin-1/2 Kondo screening. It must be stressed, how- 
ever, that after the high- 15*2 1 states are frozen out, the effec- 



tive model in the restrained Sz = ±1/2 subspace does not 
correspond to an isotropic spin-1/2 model, but rather to an 
exchange anisotropic spin-1/2 model with 



Ji_ « 2J, Jz 
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(10) 



The mapping on an anisotropic Kondo model stems from the 
spin-3/2 operators which are given in the matrix notation as 
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In the Sz = ±1/2 subspace, Sz yields ^az, while Sx.y yield 
2 X ^(Jx,y, i.e. twice the spin-1/2 operators in the transverse 
directions. The Kondo temperature is thus given by the ex- 
pression for the XXZ anisotropic Kondo model: 



Tk exp 
arctan 7 



pJz 



(12) 



which can be derived from the scaling equations to second 
order^^. For reference we also note that the expression for 
J± < Jz can be obtained by analytic continuation, giving 




Tk « Wexp 
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7 = \/l - ( ^ ) , for J_L < Jz 



These expressions somewhat overestimate the true Kondo 
temperature since they do not include the \J pJz factors; nu- 
merically calculated Kondo temperatures are shown in Fig.|3] 
The value of the Kondo temperature of the effective spin- 
1/2 screening process depends on D and it can exceed Tj^-* 
since a < 1 in the transverse dominated J ^/ Jz > case; in 
the infinite-D Hmit, Tk{S = 1/2) goes to 



Tk{S = 1/2) « 9 X IQ-^W. 



(14) 



This result fully agrees with Tk in the exchange-anisotropic 
Kondo model. Fig. [3] For finite D, the effective ratio J±/Jz 
may be different from 2, since there is a temperature range 
where the Sz = ±3/2 levels still affect the renormaliza- 
tion of the exchange interaction before they completely freeze 
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Figure 2: (Color online) Thermodynamic properties of the anisotropic Kondo model. We plot the impurity contribution to the magnetic 
susceptibility, Ximp, and the impurity contribution to the entropy. Simp. Full (blue) lines correspond to positive anisotropy, D > 0, dashed 
(red) lines correspond to negative anisotropy, D < 0, while black line is the isotropic case, D = 0. The arrow indicate in which direction the 
anisotropy is increasing. Note that the vertical axis in the upper panels is additionally rescaled as 1/S{S + 1). 



out. The intrinsic anisotropy of the effective spin- 1/2 Kondo 
model also plays a role in the splitting (shifting) of the Kondo 
resonance in a magnetic field: for strong magnetic field Bz 
in the z-direction, we expect a shift of ~ gfJ-sBz, and a 
shift approximately twice as large for fields in the transverse 
directions'^. Finally, we note that for larger spin values we 
have J± — 3 J (spin 5/2), J± = 4 J (spin 7/2), and in gen- 
eral 

J^ = {S+1/2)J. (15) 

For integer S > 2, the situation is the opposite. Now the 
regime of the exponential approach to the Sz = ground 
state corresponds to D > t\^\ For D < Tj^\ on the other 
hand, the low-temperature behavior is again that of the spin- 
1/2 Kondo model. This is made possible by the combined 
effect of the spin-S* Kondo screening at higher temperatures 
giving rise to the 5—1/2 (i.e. half-integer) residual impurity 
spin, and of the axial anisotropy which leads to freezing out 
of the high-spin Sz = ±(5* - 1/2), . . ., ±3/2 levels of this 
residual object to finally give rise to a residual spin-1/2 ob- 
ject. This residual spin-1/2 then finally undergoes a spin-1/2 
Kondo screening. At first sight it might even seem surpris- 
ing that the residual spin is compensated at all, given that in 
the isotropic high-spin Kondo model the residual exchange 
interaction is ferromagnetic, yet for D > Q the induced ex- 
change anisotropy (see Sec. HIH i is of the type which leads to 
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Figure 3: The Kondo temperature of the XXZ anisotropic 5 = 1/2 
Kondo model as a function of the Jx/Jz ratio. 

complete screening of the impurity spin*^. Since the residual 
{S ~ 1/2) spin is an extended object, there is no simple map- 
ping on the effective spin-1/2 Kondo model and it appears 
difficult to estimate the effective exchange constants J± and 
Jz even in the D limit. The effective bandwidth on the 
other hand is clearly given by W^s = aD, where a is some 
constant of order 1. The Kondo temperature is thus given 
by Tk{S = 1/2) w Df{D), where /(£>) turns out to be a 
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power-law function with non-integer exponent, f{D) ^ D^-^ 
for S* = 2 and f(D) - D^ '^ for 5 3. see Fig.g] For large 
spin t 
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Figure 4: (Color online) The cross-over temperature (indentified with 
the spin-1/2 Kondo temperature Tk for D < T^'') in integer-spin 
anisotropic Kondo models for the hard-axis case. 

We would like to emphasize the important fact that the orig- 
inal impurity spin S and the residual spin 5—1/2 after the 
Kondo screening belong to different classes (one is integer, 
the other half-integer). For this reason we always find differ- 
ent behavior depending on which of the D and is smaller. 
Furthermore, this explains similarities between integer S for 



D < T^"^ and half-integer S for D > T)^ 

Finally we consider the special case of integer spin 1 . For 
large D > rj?"*, the behavior is the same as for other inte- 
ger spins. For small D < the approach to the strong- 
coupling fixed point is of the spin-1/2 Kondo model type. 

However, we find that the second Kondo temperature de- 
pends exponentially on the D/T^' ratio, see Fig. HI a situa- 
tion strongly reminiscent of the two-stage Kondo screening^^ 
in the side-coupled impurity systems^2i^ii^2i^^i22. At T^^^ the 
spin is screened from S — 110 5= 1/2 by the spin-1 Kondo 
effect, then at T^^^ from 5 = 1/2 to 5 = by the spin-1/2 
Kondo effect. It is possible to fit the results with 
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^ K — K ^ 
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with coefficients a = 0.052, 7 = 2.531, c = 0.0661. This fit 
is only of phenomenological value: the situation here is dif- 
ferent from the one in the side-coupled impurity case, where 
it is possible to interpret the results in terms of screening of 
the second impurity by the quasiparticles resulting from the 
first stage of the Kondo screening and where Jeff corresponds 
to real exchange coupling. 

In Fig.|5]we plot the thermodynamic quantities for the re- 
lated S* = 1 Kondo model with XXZ exchange anisotropy^^. 
As expected, we find that J± < models behave similarly 
to Z? < 0, and J± > similarly to D > 0. In the presence 




Figure 5: (Color online) Thermodynamic properties of the S 
Kondo model with XXZ exchange anisotropy. 



of both symmetries, the two anisotropies may either enhance 
each other or compete. 

In Fig.|6]we plot the impurity contribution to the entropy for 
the fully anisotropic problem with the DS^ and E{S'^ — Sy) 
terms (impurity contribution to the magnetic susceptibility 
cannot be easily computed since Sz is no longer a good quan- 
tum number). We consider both the case when E is large 
[subfigure a) with E — 0.2|_D|; such ratio is found, for ex- 
ample, for integer spin-2 Fe impurities on CuN/Cu surfaces 
where D = — 1.5meV and E — 0.3 me V, and for half- 
integer spin-5/2 Mn impurities on the same surface'"] and 
for the case where E is much lower than \D\ [subfigure b) 
with E = O.OljZ?!; no such magnetic adatom/surface system 
had been identified so far in the negative-D systems, but E 
appears to be much smaller than D in the positive-D Co on 
CuN/Cu system]. 

The ground state for E ^ is non-degenerate for all pa- 
rameters. For D > 0, the 5^ — Sy operator is an iiTelevant 
perturbation and the thermodynamic behavior of the system 
is hardly affected. For D < 0, however, the perturbation is 
relevant and drives the system to a different, non-magnetic 
ground state for all S. For S = 1, the degeneracy is lifted 
on the scale of E if E > T^°^; for E < T^°^ there appears 
to be Kondo-like screening of the residual spin-1/2 with ex- 
ponential reduction of the second Kondo temperature, similar 
to what happens in the D > case. Thus the i?-term can 
induce two-stage Kondo behavior even in the D < case. 
For half-spin 5 = 3/2, the degeneracy is again lifted on the 

scale of if £; < T^"\ If > previously known 

Kondo effect with pseudo-spin 1/2 occurs^^; the Kondo tem- 
perature depends on parameters in a non-trivial way. For inte- 
ger spin 5 = 2, we find degeneracy lifting on the scale of E if 
E > and a pseudo-spin 1/2 Kondo effect if E < . 
The effective bandwidth is now given by Wcfr = aE, where 
a is some constant of order 1, so the Kondo temperature is 
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given by Tk ~ Ef{D, E) where f{D,E) has power-law be- 
havior as a function of D with a D/E (and spin S) dependent 
exponent, see Fig.|7] 

We conclude that the easy-axis systems with transverse 
anisotropy behave rather similarly to hard-axis systems; in 
the conditions for the emergence of the pseudo-spin S" = 1/2 
Kondo effect the quantity E takes the place of D. 



electron spectrum with some intermediate phase shift with ad- 
ditional splitting due to residual exchange interaction; as |Z)| 
is swept from to +oo, the (spin-dependent) phase shift is 
reduced from tt/2 to 0, while the nature of the residual im- 
purity spin changes from isotropic spin- 1/2 doublet to the 
anisotropic 5*2 = ±1 magnetic doublet, and the residual ex- 
change constants increase from to some finite value of the 
order of the bare exchange coupling J, see also App.lAl 



B. Properties of the D < systems 

For any S > I and _D < 0, £' = 0, the ground state is 
twofold degenerate. The fixed-point spectra depend on the 
value of D, thus the 5^ term is a marginal operator. At low 
temperatures, these systems have a Curie-like magnetic re- 
sponse with fractional spin 



Ximp(r^0) = (.9/iB)'^. 



(18) 



This is reminiscent of the fractional-spin non-Fermi-liquid 
fixed points in the pseudo-gap Kondo model with p (a;) ~ |cj|'' 
density of states for spin 1 and r > o^'MLSSiSl and in the re- 
lated power-law Kondo model in the case of ferromagnetic 
coupling and r < C^*. There are nevertheless some notable 
differences. The fixed points in the pseudo-gap Kondo model 
are found to be unstable with respect to the particle-hole sym- 
metry breaking. The fixed points we find are, however, sta- 
ble with respect to particle-hole symmetry breaking: both 
anisotropy, 5f , and potential scattering, Uf^, are marginal op- 
erators. Furthermore, the fixed point in the anisotropic Kondo 
model has entropy In 2 (i.e. impurity behaves as a two-level 
system), while the fixed point in the pseudo-gap Kondo model 
has entropy (1 + 2r) In 2*^"'*^'. 

In Fig. [8] we plot the Curie constant as a function of the 
D /T^^ ratio. The transition from the isotropic limit 

C^{S^ 1/2){S + l/2)/3 = - l/4)/3, (19) 

to the saturated easy-axis anisotropic behavior, 

C = 5*2, (20) 

is approximately logarithmic, i.e. very slow. This reflects the 
underlying underscreened Kondo effect and the partial screen- 
ing of the impurity moment which also feature similarly slow 
logarithmic dependence of the magnetic susceptibility on tem- 
perature and magnetic field. 

In Fig.|9]we plot the energy levels of the spin-1 negative- 
ly model as a function of the anisotropy D. For small \D\, 
the levels approach the free electron spectrum with 7r/2 phase 
shift (which is due to the Kondo effect) with additional two- 
fold degeneracy of all levels due to the decoupled spin- 1/2 
residual impurity spin. For large \D\, the levels tend towards 
a spectrum with states clustered near the energies of the free 
electron spectrum with zero phase shift (no Kondo effect) with 
additional splitting due to residual anisotropic exchange cou- 
pling that will be studied in App. lAl For intermediate \D\, 
the spectrum is more complex. It can be characterized as free 



C. Dynamic properties 

To characterize dynamic properties of the anisotropic 
models, we calculate the T-matrix and the dynamical spin 
susceptibilit y ^^-^^-^^'^^ . The T-matrix for the Kondo model can 
be determined by computing the Green's function^^ **"^'"^" 



(21) 



where — Sx i iSy are the impurity spin operators and 
/q^ creates an electron with spin cr on the first site of the hop- 
ping Hamiltonian''*'^"*. Assuming constant exchange constant 
^' = J, the T-matrix is then given by 



T{uj) J^G{u). 



(22) 



This quantity contains information on both elastic and inelas- 
tic scattering rate (cross section) ^^'^"'^^ : 



(23) 



where ctq is the scattering rate in the case of unitary 
scattering These scattering rates can, in turn, be related 
to the amplitude of the impurity -related spectral features in 
scanning tunneling spectroscopy experiments?2i2L22i21, thus in 
the following we will call the quantity o"totai(w) the "conduc- 
tance" and we will express it in the units of ctq- We also note 
that the zero-temperature scattering rate in a Fermi-liquid sys- 
tem is 



ctotaKw ^ 0) = aci{uj ^ 0) = CTo sin^((5), 



(24) 



where 6 is the quasiparticle scattering phase shift, while 

CTineKt^ ^ 0) = 0. 

The imaginary part of the dynamical spin susceptibility 
Saa{^j->) (with a = x,y, z) is defined as 

5aaH = --Im((5„;5„))^. (25) 

TT 

The dynamical spin susceptibility is in principle observable 
in tunneling-current-noise spectroscopy using spin-polarized 
STM^'*-^^. Our results thus provide information on how the 
Kondo effect in the presence of magnetic anisotropy modi- 
fies the current noise. In addition, the dynamical spin suscep- 
tibility contains information on the differential cross section 
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Figure 6: (Color online) Impurity contribution to the entropy of the anisotropic Kondo model with DS^ + E{S'^ — Sy) terms with constant 
ratio (a) E/\D\ = 0.2 and (b) E/\D\ = 0.01. We plot D = IQ^^ , 10"^, . . ., 10"^; the arrow indicate in which direction the anisotropy is 
increasins. 




10 10 
DAV 

Figure 7: (Color online) The cross-over temperature (indentified with 
the spin-1/2 Kondo temperature Tk for D < t]"') in the 5 = 2 
anisotropic Kondo model with transverse anisotropy. 



a{uj,A), where A is the energy exchange^^; the differential 
cross sections is also possibly measurable'^*. 

In Fig.[TO]we plot the dynamic properties of the anisotropic 
Kondo model with the DS^ term (E = 0). We discuss first 
the conductance. In the isotropic D = case, the conduc- 
tance always rises to the unitary limit, albeit the approach to 
this limit is slow (logarithmic)-'4i^Sdi^. Unitary scattering can 
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IDI/T, 



(0) 



Figure 8: (Color online) Curie constant as a function of the \D\/Tjp 
ratio for S = 1,3/2,2 Kondo models with easy-axis anisotropy, 
L* < 0. 



be explained by the well known fact that in the isotropic case a 
single conduction channel screens precisely 1/2 unit of the im- 
purity spin and, in doing so, the low-energy conduction band 
electrons gather a 7r/2 scattering phase shift. 

In the hard-axis D > case, the conductance at low tem- 
peratures is zero for aU integer values of the impurity spin, 
while it is unitary for all half-integer spins. For 5* = 1 and 
D < this can be explained by the two-stage Kondo 

screening: the first S — 1 screening stage leads to increased 
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Figure 9: (Color online) Energy levels at iteration A'^ = 80 (T ~ 
10~^®W) as a function of the anisotropy D in the spin-1 Kondo 
model with easy-axis anisotropy, D <Q. 



conductance as in the isotropic case, however the conduc- 

(2) 

tance drops to zero at a lower temperature scale when 
the residual spin-1/2 is compensated in the second screen- 
ing stage. This is similar to the screening of a 5 = 1 im- 
purity by two conduction channels with unequal exchange 
constants^2iS, but here the two-stage screening occurs with 
a single channel. Non-monotonous energy-dependence of the 
spectral function is also found in the case of side-coupled im- 
purities, but there both screening stages are of the S ^ \/2 
jj-jjj^50j5ij52j53 Pqj. j^j-gg ]j y t\^\ the anisotropy makes the 
impurity non-magnetic, there is no exchange scattering at en- 
ergy scales below D and the low-temperature conductance 
drops to zero. Note that in both cases the stable fixed point 
is the same, the ratio D/T^-' merely determines by which 
mechanism the 5* = 1 Kondo screening is interrupted by the 
anisotropy; the transition between the two regimes is smooth. 

For higher integer spins (represented in Fig. [TO] by the 
5 = 2 case), the situation is in some respect similar. For 
D < the conductance at first increases during the initial 
spin-S* Kondo screening, there is even a slight bump above the 
result for the isotropic case on the energy scale of D when the 
effective spin doublet is formed and there is additional reso- 
nantly enhanced scattering. Spin-S* Kondo screening yields a 
half-integer effective impurity spin, thus a spin-1/2 doublet 
on the energy scale D. The spin-1/2 Kondo screening then 
leads to a decrease in conductance to zero. As previously dis- 
cussed, the main difference from the S — \ case is that there 
is no exponential reduction of the energy scale for D < T^°\ 

For half-integer spins (represented in Fig.[TO]by the S = 
3/2 case), the conductance at low temperatures is unitary, as 

in the isotropic case. ¥oxD > T^^^this may be explained by 
the fact that the anisotropy leads in this case to a low-laying 
doublet formed by the Sz = ±1/2 magnetic doublet which 
undergoes Kondo screening like in the conventional spin-1/2 
Kondo model. The doublet is formed before the spin-S" Kondo 
screening commences and the mapping to the spin-1/2 Kondo 
model is a good approximation: the Kondo resonance then has 



approximately Lorentzian form. It may be noted that there is 
again an additional peak in (t{uj) atuj^D due to resonantly 
enhanced scattering by the Sz Sz i I processes. For D < 
the conductance curve first follows the logarithmically 
increasing isotropic spin-5 curve, then at D the approach to 
the unitary limit becomes faster. 

We emphasize the marked difference between the integer 
spin S — 2 and the half-integer spin 5 = 3/2 regarding the 
energy scale where the limiting value of (Ttotai (0 viz. 1) is 
approached: for 5 = 2, the characteristic frequency decreases 
faster than linearly with D for small D, while for 5 = 3/2 
it tends to increase very slowly with D for large D. In both 
cases this behavior is due to the underlying effective 5=1/2 
Kondo physics; for 5 = 2 it implies that the effective ex- 
change constant becomes smaller as D is reduced, while for 
S — 3/2 this is due to the saturation of the effective J±/Jz 
ratio for large D as already discussed. 

We now consider the easy-axis D < case. For all 5, 
the conductance is found to saturate at some D-dependent fi- 
nite value here. This reflects the presence of a residual un- 
compensated impurity spin which induces scattering at low 
temperatures which is neither unitary nor zero (these are the 
only two possibilities compatible with a Fermi-liquid system 
in the presence of the particle-hole symmetry which imposes 
the restriction 6 = —6 mod tt on the quasiparticle scattering 
phase shift S, thus d = or S = tt). This implies that the sta- 
ble fixed point may be classified as a non-Fermi liquid, which 
is confirmed by the non-vanishing inelastic scattering rate at 
the Fermi level, see Fig.[TT] In the limit \D\ 0, this non- 
Fermi-liquid fixed point evolves continuously into a singular 
Fermi-liquid fixed point of the isotropic underscreened Kondo 
model. Since the approach to this limit is slow (in particu- 
lar for 5 > 3/2), even very small negative axial anisotropy 
will lead to significantly reduced conductance at low temper- 
atures. In the following we show that this behavior is modified 
by non-zero E (see also Ref.[28h. Nevertheless, if i? <C \D\, 
there might exist a temperature range where the non-Fermi- 
liquid behaviour is observable. 

We now consider the imaginary part of the dynamic spin 
susceptibility Saa{^)- We recall that these quantities are 
given at zero temperature by 



1 1 



N 



\{i\Sc,\jGs)\^S{uj-{E,~EGs)) 

(26) 

where the index i runs over the excited states and the index 
Jgs over the degeneracy of the ground state, while is the 
zero-temperature partition function (equal to the degeneracy 
of the ground state). In the absence of the coupling to the 
conduction band, the transverse susceptibility 5a;a;(w) features 
delta peak(s) at 



uj = \DS^z - D{Sz - \f\ = \D{2Sz - 1)1, 



(27) 



since Sx couples neighboring levels. For 5=1, there is thus 
a peak at \D\, for 5 = 3/2 a peak at 2\D\ and for 5 = 2, 
depending on the sign of D, a peak at D (for Z) > 0) or at 
3|Z3| (for D < 0). The longitudinal susceptibility Szz{^) of a 
decoupled impurity is zero. 
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Figure 10: (Color online) Dynamic properties of the anisotropic Kondo model. We plot the spectral function of the T-matrix ("conductance") 



and the longitudinal and transversal dynamical suspectibility Szz{u}) and 5xx((jj). The curves correspond to \D\ = 10 



styles are the same as in Fig.|2] Note that for 5* ; 
scale outside the displayed range. 



. , 10 ; the curve 



2 and _D = 10 the susceptibility curves do not diverge; they start decreasing at an energy 



The S — 1 case has some special features, so we focus first 
on the regular cases, S > 3/2. The longitudinal susceptibility 
5zz((jj) always has a single peak. In the hard-axis D > case 
we find that the peak occurs at the energy scale of the effec- 
tive S = 1/2 Kondo effect (when it occurs, i.e. D > T^"* 

for 5 = 3/2 and D < t[°' for 5 = 2, respectively) or at 
~ D in the regime with no S — 1/2 Kondo effect. In both 
cases this corresponds to the energy scale where the conduc- 
tance approaches its limiting value, as discussed above. This 
result is expected, since in the presence of the Kondo effect 
the largest magnetic fluctuations always occur on the energy 
scale of TV, which is a reflection of the anomalously strong 
spin-flip scattering of the conduction band electrons off the 
impurity; in the absence of the Kondo effect, magnetic fluc- 
tuations occur on the scale of the local magnetic excitation 
energy, in this case the anisotropy D. In the easy-axis D < 
case, the position of the peak in Szz{^) is always '--^ ll?]. 

For all D and all S > 3/2, the approach to the zero fre- 
quency limit is linear, while the behavior at high frequencies 
is described approximately by a power law S{uj) ~ cj^^ with 
exponents 7 that are slightly above 1^^. 



the hard-axis D > case, we find a two-peak behavior in the 
parameter regime with S* = 1/2 Kondo effect: the first peak 
corresponds to the scale of D, the second to the scale of Tk- 
In the transition regime D ~ the two peaks merge into 
a single peak which then follows the scale of anisotropy D. 
In the easy-axis D < case, the curves again have a single 
peak at ~ \D\. At this point a comment on the peak width 
is in order: the peaks at ~ \D\ are over-broadened due to the 
broadening procedure used in the NRG method, thus the very 
narrow peaks take the form of the broadening kernel'^*. 



The transverse susceptibility Sxxii-^) is more complex. In turns negative approximately at \D\ ^ Tk , i.e. the suscep 



We now finally turn to S = 1. For hard-axis D > 
anisotropy, we find for both Szzi^^) and Sxxi^^) a single peak 
centered at w - D (for D > T^) or at w - T^^^ (for 
D < tP). In (w), there is no additional peak on the scale 
of L>forL> < tP, as was the case for the spin-2 model, but 

we observe a change of slope aX. uj ^ Even more pe- 

culiar are the results for the easy-axis D < case. The lon- 
gitudinal susceptibility has a linear frequency dependence for 
1^1 ^ Tk\ but for |D| < we observe scaling with 
exponent a which depends on the anisotropy D and which 

(0) 



11 




coAV 



Figure 1 1 : (Color online) Total, elastic and inelastic scattering cross 
sections for the 5=1 Kondo model with no anisotropy, positive 
anisotropy and negative anisotropy. 



tibility becomes divergent. In the transverse susceptibility we 
find a strong deviation from linear scaling for all values of 
|D|, with divergent behavior akeady for \D\ much larger than 

T^?' ■ We emphasize that for all spin values S, the easy-axis 
D < systems exhibit non-Fermi-liquid features (such as fi- 
nite inelastic scattering at uj ~ 0), but only for 5 = 1 is the 
dynamic spin susceptibility divergent at low frequencies. 

We note that a linear frequency dependence of Saai^^) at 
low frequencies must hold in Fermi-liquid systems as man- 
dated by the Korringa-Shiba relations^^''^''^''^ FotD < Othe 
system is, however, non-Fermi liquid and the static magnetic 
susceptibility is diverging for all S, thus S{u!) is not expected 
to be linear. We view the fact that it is linear for S* > 3/2 
merely as a coincidence. 

We finish our discussion of dynamic properties by showing 
that in the presence of the transverse anisotropy term E{S'^ — 
Sy), all non-Fermi-liquid features disappear, Fig.[TT] In par- 
ticular, we observe that for all integer S the low-temperature 
conductance is zero, while for all half-integer S it is unitary, 
irrespective of the sign and size of the axial anisotropy D. Fur- 
thermore, the spin susceptibility always has linear frequency 
dependence at low frequencies. This agrees with the results 
for static properties. 



V. CONCLUSION 

Due to the presence of the magnetic anisotropy, the im- 
purity spin is always completely compensated at low enough 
temperature, even in the case of high-spin impurities coupled 
to a single continuum channel. This goes against the common 



wisdom that a single channel can screen at most 1/2 unit of 
spin, which results in a residual spin of 5 — fc/2 in the case of 
k screening channels: such a behavior will only be observable 
for isotropic systems, a situation which is rather unlikely for 
magnetic impurities adsorbed on surfaces. 

The approach to the fully screened impurity spin with de- 
creasing energy scale depends on the impurity spin S, the 
sign of axial anisotropy D, and on the hierarchy of the energy 
scales, in particular on the dichotomy D < T^'' or _D > T^^ 
in the positive-D models. For positive D (hard-axis case), 
we find spin- 1/2 Kondo screening at low temperatures in the 
weak-anisotropy case D < for integer spin (with ex- 
ponentially reduced Kondo temperature in the special case of 
spin 1) and in the strong-anisotropy case D > for half- 
integer spin; in other cases, the effective impurity moment 
drops to zero exponentially on the energy scale of D. The 
transverse anisotropy E does not play an important role for 
D > 0. For negative D (easy-axis case), however, it changes 
the behavior of the system completely. For E = 0, the impu- 
rity spin would be only partially screened, leading to residual 
spin-dependent scattering and non-Fermi-liquid behavior. In 
real systems E ^ and the behavior of _D < becomes sim- 
ilar to that of the D > systems: for S* = 1 we again find 
two-stage Kondo screening for E < for S = 3/2 there 
is effective spin-1/2 Kondo effect for E > while for 

5 = 2 we find effective spin-1/2 screening for E < tI°\ 
If E \D\, some features of the non-Fermi-liquid behavior 
might be observable in intermediate temperature respectively 
energy regimes. The zero-temperature conductance is found 
to be zero for all integer spins and unitary for all half-integer 
spins for both signs of D (assuming E ^ 0). For integer spin 
and E,D < t'^\ the temperature dependence of the con- 
ductance is expected to be non-monotonous with increased 
conductance in the Tk{S = 1/2) < T < range. The 
most likely candidate for observation of underscreening is the 
5 = 1 case due to the exponential reduction of the second 
Kondo temperature, which is a special feature of the 5=1 
model and is not found for other integer spins. 

We conclude by noting that taking the magnetic anisotropy 
effects into account is essential in interpreting scanning tun- 
neling spectra of magnetic impurities measured at low tem- 
peratures (T ^ \D\). The appropriate low-temperature ef- 
fective impurity model depends on the parameters of the orig- 
inal physical model. While most interpretations of the sur- 
face Kondo effect have been based on the Anderson impurity 
model (which maps on the isotropic Kondo model), it is more 
likely that the appropriate effective model is, in fact, some 
anisotropic Kondo model. This is especially true for impu- 
rities on decoupling layers, but also applies to impurities ad- 
sorbed directly on metallic surface where due to stronger hy- 
bridisation charge fluctuation effects also play a role. Appro- 
priate description in terms of an anisotropic model has impor- 
tant consequences on the relation between the Kondo temper- 
ature and the "bare" model parameters (effective bandwidth, 
exchange constants and J^), and on the coupling of the 
effective impurity spin with the external magnetic field. 
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Appendix A: ANALYSIS OF THE NON-FERMI-LIQUID 
FIXED POINT 

We study the properties of the non-Fermi-liquid fixed point 
in the 5=1 Kondo model with D < in the extreme 
anisotropic limit, l^?! ^ cxd (numerical results, where needed, 
are taken from a D — — lOOM^ calculation). 

We show that in the N ^ oo limit (N being the NRG iter- 
ation number), the fixed point corresponds to fermions with 
residual anisotropic exchange scattering from a Sz = ±1 
magnetic doublet. In other words, there is no decoupling of 
the residual impurity spin, as in the isotropic 5=1 Kondo 
model^^. There is furthermore no effective potential scatter- 
ing in the \D\ oo limit, since the conduction band electrons 
cannot flip the impurity spin due to the high energy barrier, 
thus there can be no quasiparticle phase shift due to Kondo 
screening. We may thus take the Hamiltonian for the uncou- 
pled conduction band chain^ 

r 

(Al) 

where p and h are second quantization operators for particle 
and hole excitations, while £'p°'' and E^^"^ are their energies in 
the units of the characteristic energy scale at the A^-th NRG 



iteration, and ^ = In A\/A/(l - A^^) « 2.41 is an additional 
scale factor in the discretization scheme used in this work'''*. 

The fixed-point Hamiltonian can be written as H' = 
Hhand + He with 

H, = j;s,5, + Jl (s,5, + SySy) , (A2) 

where Sj are the spin operators on the first site of the Wilson 
chain, 

S^^T.fL(l<a)f0,a', (A3) 

while Si are the impurity spin-1 operators. The /o opera- 
tors need to be written in terms of the single particle and hole 
operators^: 

r 

In the \D\ — + oo limit, the Hilbert space is restricted to the 
5z = ±1 states, thus the transverse exchange drops out. 
Keeping only the terms involving the particle excitations, we 
finally obtain 

He = J* ^ aoraor'Sz'^ ~ pl^lPr'.l) ■ (A5) 

r,r' 
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Energy 


Q 




Degeneracy 








±1 


2 


0.672864 


±1 


±1/2 


4 


0.820964 


±1 


±3/2 


4 


1.345728 








2 


1.493829 


±2,0[x2] 


±1 


8 


1.641929 





±2 


2 



Table I: Excitation spectrum for large negative D (here D — 
-lOOW, A = 2.5, A'' = 79). The lowest single-particle excitations 
for free electrons is E^J'l = 0.7468559. 

The fixed-point excitation spectrum computed using NRG 
is shown in Table J] The ground state has total charge Q = 
and consists of two spin states Sz — ±1. There are four 
low-lying single-particle excited states with Q = ±1: two 
degenerate levels with Sz = ±1/2 and (scaled) energy 

E,{A) = ~ ^^>oiA(^-^)/^e (A6) 

and two other degenerate levels with Sz — ±3/2 and energy 

E,iB) = El^l + y:al,A^^-'^^'^ (A7) 



In the first two levels, the impurity spin and the spin of the 
particle excitation are anti-aligned, while in the second two 
levels the impurity spin and the spin of the particle are aligned. 

Since the energy difference between Ep{A) and Ep{B) 
no longer varies with N at the fixed point and the factor 
all — Q!oi A^^"^)/^ reaches its Umiting value (equal to 0.309 
for A = 2.5) this implies that J* is a constant given by 

j: = ^A^l^IA^. (A8) 

For pj = 0.1 and D = -lOOW (TableU we obtain 

j; w 0.198M^ = 0.99J. (A9) 

As expected, we obtain J* which is essentially equal to the 
bare J in the large- |D| limit. 



ACKNOWLEDGMENTS 

We thank A. Heinrich for sharing unpublished results for 
tunneling spectra of Co impurities. We acknowledge sup- 
port by Gesellschaft fiir wissenschaftliche Datenverarbeitung 
(GWDG) in Gottingen through SFB 602 (T. P. and R. Z.) and 
project PR 298/10-1 (R.P). 



* S. Heinze, M. Bode, O. Pietzsch, A. Kubetzka, X. Nie, S. Blugel, 
and R. Wiesendanger, Science 288, 1805 (2000). 

^ C. A. R Vaz, J. A. C. Bland, and G. Lauhoff, Rep. prog. phys. 71, 
056501 (2008). 

^ R. Wiesendanger, 1. V. Shvets, Diirgler, G. Tarrach, H. J. 
Guntherodt, J. M. D. Cody, and S. Graser, Science 255, 583 
(1992). 

■* M. Farle, Rep. Prog. Phys. 61, 755 (1998). 

^ P. Gambardella, A. Dallmeyer, K. Maid, M. C. Malagoli, W. Eber- 
hardt, K. Kern, and C. Carbone, Nature 416, 301 (2002). 

* S. Rusponi, T. Cren, N. Weiss, P. Buluschek, L. Claude, and 
H. Brune, Nature Materials 2, 546 (2003). 

W. Kuch, Nature materials 2, 505 (2003). 
^ A. J. Heinrich, J. A. Gupta, C. P. Lutz, and D. M. Eigler, Science 
306, 466 (2004). 

' M. Bode, M. Heide, K. von Bergmann, P. Ferriani, S. Heinze, 
G. Bihlmayer, A. Kubetzka, O. Pietzsch, S. Bliigel, and 
R. Wiesendanger, Nature 447, 190 (2007). 

C. F Hirjibehedin, C.-Y. Lin, A. F Otte, M. Ternes, C. P Lutz, 
B. A. Jones, and A. J. Heinrich, Science 317, 1199 (2007). 
" P. Ferriani, K. von Bergmann, E. Y. Vedmedenko, S. Heinze, 
M. Bode, M. Heide, G. Bihlmayer, S. Bliigel, and R. Wiesendan- 
ger, Phys. Rev. Lett. 101, 027201 (2008). 

P. Gambardella, S. Rusponi, M. Veronese, S. S. Dhesi, C. Grazioli, 
A. Dallmeyer, I. Cabria, R. Zeller, P. H. Dederichs, K. Kern, et al.. 
Science 300, 1130 (2003). 

A. C. Hewson, The Kondo Problem to Heavy- Fermions (Cam- 
bridge University Press, Cambridge, 1993). 
R. Bulla, T. Costi, and T. Pruschke, Rev. Mod. Phys. 80, 395 
(2008). 



O. Ujsaghy, A. Zawadowski, and B. L. Gyorffy, Phys. Rev. Lett. 
76, 2378 (1996). 

O. Ujsaghy and A. Zawadowski, Phys. Rev. B 60, 10602 (1999). 
" L. Szunyogh, G. Zarand, S. Gallego, M. C. Mufioz, and B. L. 

Gyorffy, Phys. Rev. Lett. 96, 067204 (2006). 
"* A. Szilva, S. Gallego, M. C. Muiioz, B. L. Gyorffy, G. Zarand, 

and L. Szunyogh, Friedel oscillations induced surface magnetic 

anisotropy, arxiv:0805.3275 (2008). 
" C. Etz, J. Zabloudil, P. Weinberger, and E. Y. Vedmedenko, Phys. 

Rev. B 77, 184425 (2008). 

V. Madhavan, W. Chen, T. Jamneala, M. Crommie, and N. S. 
Wingreen, Science 280, 567 (1998). 
^' J. Li, W.-D. Schneider, R. Berndt, and B. Delley, Phys. Rev. Lett. 
80, 2893 (1998). 

O. Ujsaghy, J. Kroha, L. Szunyogh, and A. Zawadowski, Phys. 
Rev. Lett. 85, 2557 (2000). 

A. Schiller and S. Hershtield, Phys. Rev. B 61, 9036 (2000). 
V. Madhavan, W. Chen, T. Jamneala, M. F. Crommie, and N. S. 
Wingreen, Phys. Rev. B 64, 165412 (2001). 
A. J. Heinrich, private communication. 

H. B. Heersche, Z. de Groot, J. A. Folk, H. S. J. van der Zant, 

C. Romeike, M. R. Wegewijs, L. Zobbi, D. Barreca, E. Tondello, 
and A. Comia, Phys. Rev. Lett. 96, 206801 (2006). 

M. Jo, J. Grose, K. Baheti, M. Deshmukh, J. Sokol, E. Rumberger, 

D. Hendrickson, J. Long, H. Park, and D. Ralph, Nano Lett. 6, 
2014 (2006). 

C. Romeike, M. R. Wegewijs, W. Hofstetter, and H. Schoeller, 
Phys. Rev. Lett. 96, 196601 (2006). 
^' C. Romeike, M. R. Wegewijs, W. Hofstetter, and H. Schoeller, 
Phys. Rev. Lett. 97, 206601 (2006). 



14 



D. Roosen, M. R. Wegewijs, and W. Hofstetter, Phys. Rev. Lett. 

100, 087201 (2008). 
" L. Bogani and W. Wernsdorfer, Nature materials 7, 179 (2008). 
'2 D. M. Cragg and R Lloyd, J. Phys. C: Solid State Phys. 12, L215 

(1979). 

D. M. Cragg, P Lloyd, and P Nozieres, J. Phys. C: Solid St. Phys. 

13, 803 (1980). 

P. Mehta, N. Andrei, P. Coleman, L. Borda, and G. Zarand, Phys. 
Rev. B 72, 014430 (2005). 

W. Roller, A. C. Hewson, and D. Meyer, Phys. Rev. B 72, 0451 17 
(2005). 

K. Furuya and J. H. Lowenstein, Phys. Rev. B 25, 5935 (1982). 
N. Andrei, K. Furuya, and J. H. Lowenstein, Rev. Mod. Phys. 55, 

331 (1983). 

N. Andrei and C. Destri, Phys. Rev Lett. 52, 364 (1984). 

N. Andrei, Integrable models in condensed matter physics, cond- 

mat/9408011 (1994). 

P. D. Sacramento and P. Schlottmann, Phys. Rev B 40, 431 
(1989). 

P. D. Sacramento and P. Schlottmann, J. Phys.: Condens. Matter 

3, 9687 (1991). 



P Nozieres and A. Blandin, J. Physique 41, 193 (1980). 
P Coleman and C. Pepin, Phys. Rev B 68, 220405(R) (2003). 



43 

M. Bortz and A. Klumper, Eur. Phys. J. B 40, 25 (2004). 
'^^ D. C. Mattis, Phys. Rev Lett. 19, 1478 (1967). 

P Schlottmann, Phys. Rev. Lett. 84, 1559 (2000). 

P Schlottmann, J. Appl. Phys. 89, 7183 (2001). 
'^^ A. Posazhennikova and P Coleman, Phys. Rev. Lett. 94, 036802 

(2005). 

A. Posazhennikova, B. Bayani, and P. Coleman, Phys. Rev. B 75, 
245329 (2007). 

5° M. Vojta, R. Bulla, and W. Hofstetter, Phys. Rev B 65, 140405(R) 

(2002). 

P S. ComagUa and D. R. Grempel, Phys. Rev B 71, 075305 
(2005). 

" R. Zitko and J. BonCa, Phys. Rev B 73, 035332 (2006). 

C.-H. Chung, G. Zarand, and P Wolfle, Phys. Rev B 77, 035120 

(2008). 

R. Zitko and J. Bonca, Phys. Rev B 74, 045312 (2006). 
R. Zitko and J. Bonca, Phys. Rev B 76, 241305(R) (2007). 
^'^ W. Izumida, O. Sakai, and S. Tarucha, Phys. Rev Lett. 87, 216803 
(2001). 

M. Pustilnik and L. I. Glazman, Phys. Rev Lett. 87, 216601 
(2001). 

M. Pustilnik and L. I. Glazman, Phys. Rev B 64, 045328 (2001). 

O. Sakai and W. Izumida, Physica B 328, 125 (2003). 
*° M. Pustilnik, L. 1. Glazman, and W. Hofstetter, Phys. Rev. B 68, 

161303(R) (2003). 
" W. Hofstetter and G. Zarand, Phys. Rev. B 69, 235301 (2004). 

J. Schmid, J. Weis, K. Eberl, and K. v Klitzing, Phys. Rev. Lett. 

84, 5824 (2000). 

S. Sasaki, S. de Franceschi, J. M. Elzerman, W. G. van der Wiel, 
M. Eto, S. Tarucha, and L. P. Kouwenhoven, Nature 405, 764 
(2000). 

^ W. G. van der Wiel, S. De Franceschi, J. M. Elzerman, S. Tarucha, 
L. P. Kouwenhoven, J. Motohisa, F. Nakajima, and T. Fukui, Phys. 



Rev Lett. 88, 126803 (2002). 

A. Kogan, G. Granger, M. A. Kastner, D. Goldhaber-Gordon, and 
H. Shtrikman, Phys. Rev B 67, 113309 (2003). 
A. Fuhrer, T. Ihn, K. Ensslin, W. Wegscheider, and M. Bichler, 
Phys. Rev. Let. 91, 206802 (2003). 

G. Granger, M. A. Kastner, I. Radu, M. P. Hanson, and A. C. 
Gossard, Phys. Rev. B 72, 165309 (2005). 

R. M. Konik, H. Saleur, and A. W. W. Ludwig, Phys. Rev. B 66, 

075105 (2002). 

A. Schiller and L. De Leo, Phys. Rev. B 77, 075114 (2008). 
™ P W. Anderson, J. Phys. C: Solid St. Phys. 3, 2436 (1970). 
^' A. M. Tsvelick and P B. Wiegmann, Adv Phys. 32, 453 (1983). 

T. A. Costi and C. Kieffer, Phys. Rev. Lett. 76, 1683 (1996). 
" T. A. Costi, Phys. Rev. Lett. 80, 1038 (1998). 
■'^ K. G. Wilson, Rev Mod. Phys. 47, 773 (1975). 

H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. 
Rev B 21, 1003 (1980). 



D. Kaihe and B.-H. Zhao, Phys. Lett. A 340, 337 (2005). 



76 

N. Andrei and J. H. Lowenstein, Phys. Rev. Lett. 46, 356 (1981). 

C. Jayaprakash, H. R. Krishna-murthy, and J. W. Wilkins, Phys. 
Rev. Lett. 47, 737 (1981). 

™ R. Peters and T. Pruschke, New Journal of Physics 8, 127 (2006). 
^° C. Gonzalez-Buxton and K. Ingersent, Phys. Rev B 57, 14254 
(1998). 

^' S. Horens and M. Vojta, Phys. Rev. B 72, 115117 (2005). 

D. Withoff and E. Fradkm, Phys. Rev Lett. 64, 1835 (1990). 

K. Chen and C. Jayaprakash, J. Phys.: Condens. Matter 7, L491 
(1995). 

M. Vojta and R. Bulla, Eur. Phys. J. B 28, 283 (2002). 
T. A. Costi, A. C. Hewson, and V. Zlatic, J. Phys.: Condens. Mat- 
ter 6, 2519(1994). 

W. Hofstetter, Phys. Rev Lett. 85, 1508 (2000). 

R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev B 74, 2451 14 

(2006) . 

A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 076402 

(2007) . 

T. A. Costi, Phys. Rev. Lett. 85, 1504 (2000). 
* G. Zarand, L. Borda, J. von Delft, and N. Andrei, Phys. Rev. Lett. 

93, 107204 (2004). 
" L. Borda, L. Fritz, N. Andrei, and G. Zarand, Phys. Rev B 75, 

235112 (2007). 

M. Vojta, R. Zitzler, R. Bulla, and T. Pruschke, Phys. Rev. B 66, 
134527 (2002). 

P S. Cornagha and C. A. Balseko, Phys. Rev B 67, 205420 

(2003). 

A. V. Balatsky, Y. Manassen, and R. Salem, Phys. Rev. B 66, 
195416 (2002). 

Z. Nussinov, M. F. Crommie, and A. V. Balatsky, Phys. Rev. B 68, 
085402 (2003). 

M. Garst, P. Wolfle, L. Borda, J. von Delft, and L. Glazman, Phys. 
Rev. B 72, 205125 (2005). 



97 



H. Shiba, Prog. Theor. Phys. 54, 967 (1975). 



W. Hofstetter and S. Kehrein, Phys. Rev B 63, 140402(R) (2001). 
V. L. Campo and L. N. Oliveira, Phys. Rev. B 72, 104432 (2005). 



